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I. INTRODUCTION 

Since its discovery in 1968 [lj . the Lax pair formula- 
tion has played an invaluable role in studying integra- 
bility of various systems. Although first formulated for 
systems with infinite degrees of freedom, the formalism 
can be also used for, and provides an elegant description 
of, special finite dimensional systems with symmetries. 
Examples of such (completely integrable) systems admit- 
ting a Lax pair formulation include the Kepler problem, 
the Euler, Lagrange and Kowalevski tops, the Neumann 
model, or Toda lattice; we refer the reader to the mono- 
graph and to references therein. In what follows we 
concentrate on finite dimensional systems. Namely, we 
shall discuss a geometrization of the Lax pair matrices 
for motion in curved spacetime. 

A standard dynamical system is described on a phase 
space P equipped with the symplectic 2-form VL and with 
the corresponding Poisson brackets {,}. The dynamics 
is encoded in a Hamiltonian H through the evolution 
equation for an arbitrary scalar observable F, 



F = {F,H}. 



(1.1) 
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The Lax pair method consists of finding two phase- 
space valued matrices L and M, such that the equations 
of motion imply the Lax pair equation 



L = [L,M] 



(1.2) 



The stronger formulation requires that Eq. (|1.2[) implies 
the equations of motion, in which case the Lax pair for- 
mulation can be used as a starting point of the descrip- 
tion of the dynamical system. 

However, even without this latter stronger condition, 
the Lax pair matrices satisfying (|1.2j) play an important 
role in the study of integrability since they allow a simple 
construction of constants of motion. Indeed, the solu- 
tion of (T2J) is of the form L(t) = G(t)L(0)G- 1 (i), where 
the evolution matrix G(t) is determined by the equation 
G = — MG. Therefore, if I(L) is a function of L invariant 
under conjugation L — > GLG -1 , then I(L(t)) is a constant 
of motion. All such invariants can be generated from the 
traces of various matrix-powers of L: 



tr(L J 



(1.3) 



The particular Lax pair may not yield all the constants 
of motion. However, in such a case it is often possible to 
upgrade the initial Lax pair so that the upgraded one al- 
ready yields all the conserved observables of the dynami- 
cal system. Since the dimensionality of the Lax matrices 
is not fixed and the Lax pair equation is linear, two Lax 
pairs can be easily combined by their direct sum. An- 
other useful method of producing a parametric class of 
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Lax pairs is to introduce so-called spectral parameters, 
see, e.g., @. 

Unfortunately, in general there is no constructive pro- 
cedure to find a Lax pair for the given problem or even 
to determine whether the Lax pair (in its stronger for- 
mulation) exists. Moreover, the solution is in no sense 
unique and even the dimensionality of the matrices may 
vary. However, when the Lax pair exists, it can be a very 
powerful tool for dealing with the conserved quantities. 

In our paper, we focus on the construction of the Lax 
pair for motion on a curved manifold. In this case the 
phase space P is given as a cotangent bundle of a configu- 
ration space M. There are thus two important features of 
such a theory which single it out among generic dynam- 
ical systems: (i) It has preferred splitting of 2D phase 
coordinates into two sets — D spacetime coordinates x a 
and D momenta p a . (ii) The configuration space has 
usually an additional structure on it, the metric g a b- 

The prominent example of such dynamical systems is 
motion of particles and light in curved spacetime which 
plays an important role in General Relativity and its 
generalizations. We will concentrate mainly on this sys- 
tem. However, the dynamical systems with the cotangent 
bundle structure include also all non-relativistic systems 
which start with the Lagrangian description on the con- 
figuration space. The metric on the configuration space 
in such cases emerges from the kinetic part of the Hamil- 
tonian. 

Our aim is to show that for these dynamical systems it 
is possible to formulate a covariant analogue of the Lax 
equation. This covariant formulation allows us to em- 
ploy the description of explicit and hidden symmetries of 
curved manifolds encoded in structures as Killing vectors, 
Killing tensors or Killing- Yano tensors. In particular, for 
motion of particles in a curved spacetime we establish 
a relationship between conserved quantities, connected 
with the Lax pair, and integrals of motion connected with 
hidden symmetries, generated by Killing- Yano tensors. 

To achieve this program, we first (in Sec. [TTJ) analyze 
geometrical structures on the phase space with cotan- 
gent bundle structure, especially those induced from a 
configuration-space covariant derivative. We define a 
covariant derivative acting on phase-space fields with 
configuration-space tensor indices. 

With this geometrical background, in Sec. IIIII we for- 
mulate the covariant Lax equation 



7t La > = 



(1.4) 



for a phase-space dependent tensor field L a b (x,p) and call 
this object a Lax tensor. Here, tt is the time derivative 
(the derivative along the Hamiltonian flow) defined in 
terms of the covariant derivative mentioned above. 

By construction, any scalar invariant built covariantly 
from the Lax tensor (or from a set of Lax tensors) is pre- 
served along phase space trajectories and defines thus a 
constant of motion. Moreover, all scalar invariants en- 
coded in one Lax tensor L% can be generated from the 



traces of powers of L%, 

tr(U) =L\L\...L d a 



(1.5) 



total of j tensor: 



which is an obvious geometric analogue of Eq. (|1 .3[) . 

There is a direct translation of the covariant Lax equa- 
tion (|1.4p into the standard Lax pair equation (|1.2|) . 
Namely, the components of Lax tensor L a \ ) form the Lax 
matrix L, while the other matrix M is given through the 
(Christoffel) connection symbols as follows: 



M 



(1.6) 



An analogous geometric construction using Clifford ma- 
trices will be also introduced. Similar ideas in special 
cases have been already studied in works ■ 

The paper is organized as follows. In the next sec- 
tion we review the phase space formalism for the phase 
space built as a cotangent bundle of a configuration man- 
ifold. In particular, we define the covariant phase space 
derivative. This is further elaborated in the Appendix. In 
Sec. lIIII we recapitulate the Lax pair formalism and derive 
its geometric covariant counterpart encoded in Eq. (|1.4[) . 
Sec. IIVI illuminates the previous discussion by studying 
two 'trivial' examples of Lax tensors for geodesic mo- 
tion in generic spacetimes. Sec. |V] is devoted to highly 
non-trivial examples of Lax tensors for geodesic motion 
in special spacetimes admitting a hidden symmetry of 
Killing- Yano tensors. The motion of a charged particle 
in weakly charged Kerr-NUT spacetimes is discussed in 
Sec. ED We conclude in Sec. EED 



II. COVARIANT DERIVATIVE ALONG A 
PHASE SPACE TRAJECTORY 

Phase space 

Particle motion in a curved spacetim^B M can be de- 
scribed in the language of Hamiltonian mechanics. The 
phase space P is the cotangent bundle T*M equipped 
with the standard symplectic structure f2 and Poisson 
brackets. A phase space point can be written as [x,p] 
with position x € M and momentum p € T* X M . Any 
spacetime coordinates x a together with the correspond- 
ing components p a of the momentum p form the canonical 
coordinates x a , pb in which the symplectic structure and 



As we mentioned in Introduction, although we aim our discus- 
sion mainly to the relativistic context, the formalism introduced 
here does not depend on the signature of the metric and it can 
be used also in the case of standard non-relativistic mechanics 
with a curved configuration space. We use small Latin letters 
for spacetime (configuration space) tensor indices and we drop 
indices when it does not cause confusion. 
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Poisson brackets read 



n = dx a Adp a , (2.1) 

IF r\ = ?L™-?L^- (22) 
1 ' r dx a d Pa d Pa dx a ■ v ' ' 

Given a Hamiltonian H, the time derivative of any 
observable F is 



F = {F, H} , 



(2.3) 



which can be interpreted also as the derivative along the 
Hamiltonian flow Xj^ = fl^ 1 ■ dH , 



F = X H -dF 



(2.4) 



See Eq. (|A22[) and Fig. Q] in Appendix for further discus- 
sion. 

The covariant derivative V^- along the direction 
X = [u, f] acting on phase-space fields with spacetime 
indices is defined by the following rules: 

(i) For a field depending only on the spacetime posi- 
tion, A b y(x,p) = a b '"(x), the derivative reduces to 
the standard covariant derivative along u: 

(ii) For the momentum field p a , the derivative gives the 
momentum part of X: 

Vjf Pa = fa ■ 



More details on the conventions used for the symplectic 
structure and its inverse can be found in the Appendix. 



Covariant phase space derivative 

The derivative (j2.4[) along Xjj is defined only for scalar 
phase-space observables. It turns out fruitful to gen- 
eralize it to more general observables, namely to fields 
on phase space with spacetime indices. Such fields ap- 
pear naturally as a combination of spacetime tensors con- 
tracted with momenta, e.g., 



A a b ;;;(x, P ) = a a b -;(x) + b a b ;;; k (x) Pk + ^ 



M 



(x) PkPl 



(2.5) 

One can also consider fields of more complicated analytic 
form, e.g., A a (x,p) = (g kl (x)p k pi)~ 5 p a . 

To define a derivative of tensor fields one needs an ad- 
ditional structure. For spacetime fields such a structure 
is a covariant derivative. Assuming the covariant deriva- 
tive V on the spacetime we lift this derivative to act on 
phase-space fields with spacetime indices. We call the 
resulting operation a covariant phase space derivative. 

This derivative is defined along a general phase-space 
direction X 6 TP. The direction can be represented us- 
ing a configuration direction u £ TM and a momentum 
direction / 6 T*A/. The definition of this splitting re- 
quires the covariant derivative V. Namely, if X is tangent 
to a phase-space trajectory [x(t),p(t)], its configuration 
part u is tangent to the configuration trajectory x(t) and 
the momentum part / is the covariant derivative of p(t) 
along x(t). That is, we have 



fa 



dt 



Pa = Pa 



Pk 



(2.6) 



(iii) The derivative Vj satisfies all standard rules for 
the derivative (i.e., linearity, the Leibniz product 
rule, and the chain rule). 

These rules reflect the splitting of the phase-space di- 
rection into configuration and momentum parts. This 
can be also encoded using the partial derivatives!! and 
introduced in Appendix: 

^:::-^;;; + f n ^-A t ;;;. (2.7) 

For a phase-space trajectory with tangent field X we nat- 
urally write 



—At" = V X AI 



dt 



(2.8) 



and call ^ the covariant derivative along the phase space 
trajectory or just the covariant time derivative. The con- 
figuration and momentum parts of the Hamiltonian flow 



Xjj are 



d Pa > u 



dx 



cf. Eq. (|A7|) . which gives 



V n H d 



V A a..._dH V . „ 

dt b - 8p n Ox b - dx d Pn b - 



(2.9) 



(2.10) 



This is a natural generalization of Eq. (|2.3|) to the case 
of tensor- valued observables. 

The introduced phase-space covariant derivative can 
be expressed in coordinates, 



W x At;;; = At:::+ u n T^At; 



n-pfe «... 



(2.11) 



2 This can bo an arbitrary covariant derivative, not necessarily the 
metric one. For simplicity, in the body of the paper we assume 
vanishing torsion. Expressions with torsion can be found in the 
Appendix. In the case of geodesic motion we choose the metric 
covariant derivative. 



3 In short, the partial derivative is the derivative in a momen- 
tum direction with x fixed and ^ is the derivative in a config- 
uration direction with p parallel-transported. For more details, 
see flAH and <A2t . 
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where A^- - is the derivative of components of A along X 
direction, 



ia. 
A b.. 



dx T ' 



dp n 



(2-12) 



Note that the coordinate time derivative p a and the co- 
variant derivative f a are related by Eq. (|2.6I) . 

For more details on the derivatives V^- , ^ , and the 
corresponding coordinate expressions we refer the reader 
to the Appendix. 



Derivative of Clifford fields 

The spacetime metric derivative can also be lifted to a 
derivative acting on phase-space fields with Dirac spinor 
indices. Namely, we are interested in Clifford objects, i.e., 
operators acting on Dirac spinors. These are generated 
by the abstract gamma matrices 7°, obeying^ 



a b , b a n ab-n 

77+77 = 2g 1 



(2-13) 



Thanks to this rule, any Clifford object ([) can be 
represented by an inhomogeneous antisymmetric form 

w = Er Tuj > 



^ = £^a,.. Q .7 ai -^ 



(2-14) 



orthonormal frame ea G TM and the spinor frame Ey 
in such a way that the components!! 7° of the gamma 
matrices are constants. The covariant derivative of the 
Dirac spinor $ expressed in the spinor frame then reads 



(2-18) 



with the connection coefficients S a uniquely determined 
in terms of the Ricci coefficients T£l = (V a e|) e™ by the 
standard relation 



1 



(2.19) 



With these definitions the covariant derivative of the 
Clifford field A(x,p) on the phase space is 



V X A = A+[u n £ n , A] , 



(2.20) 



where A is just an ordinary derivative along X of com- 
ponents A*^. 



III. COVARIANT LAX EQUATION 



Lax pair 



Here, r co are homogeneous rank-r antisymmetric forms 
and 



Oi...a r _ [ai ,..^a f 



(2.15) 



The covariant derivative on the Dirac bundle is induced 
from the spacetime metric derivative by the condition^ 



V„7 Q = 



(2.16) 



Clearly, we can lift this derivative to act on the Clif- 
ford valued fields on the phase space in a similar way as 
we did for spacetime-tensor valued fields. Namely, for 
A(x,p) = \ ai ...a r (x,P) -f ai - a *-( x ) we simply get 



V X A= (V x A ai ... ar ) 7 Q1 - a " ■ 



(2.17) 



To write down this covariant derivative in components, 
in addition to coordinates x a , one has to introduce an 



4 In expressions with Clifford objects and Dirac spinors the Clif- 
ford multiplication is assumed. In components, it reduces to the 
standard matrix multiplication. 

5 To define the covariant derivative on Dirac spinors uniquely, the 
condition 112.1611 must be supplemented by some further condi- 
tions reflecting irrcducibility and reality properties of 7™. How- 
ever, thanks to Eq. 112. 161 ) and rule I I2.13H . only the covariant 
derivative which annihilates the metric can be lifted to the Dirac 
bundle. In this context we always assume vanishing torsion, so 
the derivative on tangent bundle must be the metric derivative. 



The Lax pair provides a useful tool for generating con- 
served quantities. The phase-space valued matrices L and 
M form the Lax pair if they satisfy the Lax pair equa- 
tion Q 



[L, Ml 



(3.1) 



Here, the dot is understood as the ordinary time deriva- 
tive of each component of the matrix L. It follows that 
any scalar invariant formed from the matrix L is a con- 
served quantity, cf. Eq. (|1.3p . 

It is customary to require that the Lax pair satisfies ad- 
ditional properties, especially, that (i) Eq. (|3.1[) is equiv- 
alent to the equations of motion, and that (ii) the in- 
variants of L generate the maximal number of conserved 
quantities of the system. However, we will study Lax 
pairs without requiring these additional conditions. This 
is justified by realizing that the Lax-pair equation (|3.ip is 
linear in both L and M. One can thus obtain a 'more so- 
phisticated' Lax pair as a direct sum of smaller matrices, 
each of which satisfy Eq. (]3 . 1 [) . and impose additional 
conditions only at the end, on the resulting pair. 



6 Components with respect to the frame e„ will be denoted with 
hatted indices. We will mostly skip the spinor indices (capital 
Greek letters), i.e., instead of 7 a ** we write just 7". The matrix 
multiplication between Clifford objects and spinors is assumed. 
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Lax tensor 



Clifford Lax tensor 



Wc shall now formulate an alternative covariant de- 
scription of the Lax equation and clarify its relation to 
the standard Lax pair formulation. Using the above def- 
inition of the covariant phase space derivative, we define 
the Lax tensor to be a covariantly conserved tensor field 
L a b (x,p), obeying the covariant Lax tensor equation 



V 

di 



Ll = 



(3.2) 



Any covariantly conserved_| antisymmetric form A a f,.. 
induces the Clifford field A = A a f,...7 ab " satisfying 



V 

di 



A = 0. 



(3.5) 



We call such A a Clifford Lax tensor. In components, 
using Eq. ([2~2T)|) . Eq. fl33|) implies 



A= [A, u*E n ] 



(3.6) 



Obviously, any scalar covariantly constructed^ from the 
Lax tensor or a set of Lax tensors is a constant of motion. 
In particular, this is true for invariants constructed as 
traces (|1.5p of various powers of L a . Similar to the Lax 
pair, the Lax tensor hence generates constants of motion. 

Moreover, each Lax tensor defines a Lax pair. Indeed, 
in components, while using (|2.11[) . Eq. (|3.2p implies 



L% = Ll 



%j^k nj^a j k 

L nh U 1 , nk Lj b 



(3.3) 



This equation is in form already very close to (|3.1I) . How- 
ever, the matrices forming the Lax pair must be defined 
as functions on phase space. Therefore we have to elimi- 
nate the velocity u using the first of the Hamilton equa- 
tions (|2.9p . The correspoding Lax pair matrices thus are 



[L a b ] 



M 



(3.4) 



Hence, equation (|3.2j) can be understood as a 'covariant 
generalization' of the Lax pair equation (|3.1[) ; its coordi- 
nate form gives the Lax pair in the ordinary sense. 

For a motion in curved space governed by a given 
Hamiltonian, the covariant Lax tensor satisfying p.2[) 
and the ordinary Lax pair matrices (|3.4p carry the same 
information. The Lax pair matrices are, however, co- 
ordinate dependent. Under a change of coordinates the 
matrix L transforms just by a trivial conjugation. How- 
ever, the matrix M changes in a more complicated man- 
ner since the connection coefficients are involved. 



7 The covariant derivative employed in the Lax equation 1 13.211 can 
be arbitrary. By covariant construction then we mean any op- 
eration which commutes with this derivative. Typical covariant 
operations are traces, contracted multiplication, or the determi- 
nant. If a tensor covariantly constant with respect to the cho- 
sen derivative is available, it can be used to construct the con- 
served scalar. Therefore, we typically choose the metric covari- 
ant derivative since then the metric "is available" . However, one 
could use a different derivative, for example, if the Hamiltonian 
is H = ^p a pbk ab with the "inverse mass" tensor k ab different 
from the metric g ab . 



In a way similar to previous subsection we find that ma- 
trices 



L = M = [A aS 



.7 



M 



dH 

OPn 



(3.7) 



form a Lax pair satisfying (|3.1[) . 

In the next two sections we shall give a number of ex- 
amples of (Clifford) Lax tensors for geodesic motion in 
curved spacetimc. Whereas the following section concen- 
trates on 'trivial' examples in generic spacetimes (no en- 
hanced symmetry is assumed), in Sec. |V]we discuss Lax 
tensors in special spacetimes admitting hidden symme- 
tries. The motion of a charged particle in weakly charged 
Kerr-NUT spacetimes is discussed in Sec. ED 



IV. LAX TENSORS AND GEODESIC MOTION: 
TWO TRIVIAL EXAMPLES 

Geodesic motion 

A geodesic motion with respect to the spacetime metric 
g a b is governed by the Hamiltonian 



H=^- P a9 ab Pb 

2m 



(4.1) 



The equations of motion (|2.3|) for canonical coordinates 
x a , p a are 



x a = m" V>n 



A7 



Pa 



2m dx a 



PkPi ■ 



(4.2) 



These equations are equivalent to the geodesic equation 
it n V„u a = with the velocity u a = x a and V being the 
metric covariant derivative. They are also equivalent to 
the covariant equations on the phase space: 



Pa = mg an u , 



(4.3) 



8 In the context of Clifford fields we always assume the metric 
covariant derivative since the abstract gamma matrices are then 
covariantly constant 12.161 . cf. also footnote [7] 
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Lax tensor implying geodesic motion: Example I 



As a first example, let us study the simple tensor 



t " 



9 an PnPb ■ 



(4.4) 



Thanks to Eq. (|4.3[) , it is covariantly conserved and hence 
defines a Lax tensor. 

In fact, we shall prove now that the existence of this 
Lax tensor is equivalent to the geodesic motion and hence 
it can be used as a starting point of the dynamics. In 
order to do that, we investigate the following problem: 
Assuming that the covariant Lax equation (|3.2[) along an 



unknown Hamiltonian flow Xjj is satisfied for (|4.4j) . does 
it imply that the Hamiltonian H must generate geodesic 
motion? 

Substituting (|L4|) into (pOj) . we obtain 



dt 



L a b = r Pb + P a f b = o 



(4.5) 



which for generic momentum p implies / = 0. However, 
for the Hamiltonian flow, / is given by (|2.9[) and we have 



VH 

dx 







(4.6) 



The Hamiltonian must thus be constructed only from 
momenta and covariantly constant spacetime tensors. 
In a generic curved spacetime the only covariantly con- 
stant spacetime tensors are constructed from the metric. 
Therefore, the Hamiltonian must be of the form 



H 



1 



Hp 2 



V = PaPb9 



ah 



(4.7) 



the Lax pair equation (|3. II) implies the following condi- 
tion on the Hamiltonian: 



. 



(4.10) 



Obviously, the geodesic Hamiltonian H = -^p 2 solves 
this condition. However, it is not clear if this is a unique 
solution. 



Lax tensor implying geodesic motion: Example II 

Allowing for Clifford-valued fields, it is possible to 
write an even simpler Lax tensor for the geodesic mo- 
tionfl 



A 



Pal 



(4-11) 



Similar to the Example I, this Lax tensor exists in a 
generic spacetime, it generates only conserved quantities 
which arc functions of p 2 , and the corresponding covari- 
ant Lax equation is equivalent to the geodesic equations 
of motion. 

Using (|3.7p . the corresponding Lax pair matrices are 



[A] 



r l 



P Z-iJ 



[PaT 
1 

Am 



P 



L nkll 



(4.12) 



Similarly to Example I, it is not clear to us whether these 
Lax pair matrices imply geodesic motion. Namely, the 
commutator in the Lax equation (|3.ip can be simplified 
using properties of the gamma matrices 



1 



PaT 



_P Prisma I 



(4.13) 



with h being an arbitrary function. It implies the velocity 
u a = h' (p 2 )g an p n - Since p 2 is conserved, we reproduced 
the geodesic equations of motion. 

The significance of the Lax tensor (|4.4j) lies thus in the 
fact that it implies geodesic motion. On other hand, the 
only conserved quantities which can be obtained from 
this Lax tensor are functions of p 2 (which is, of course, 
the only conserved quantity for generic geodesic motion). 

Using (|3.4[) . the Lax tensor (|4.4[) defines the following 
Lax pair matrices: 



L% = 9 an P nPb 



1 



b = —Pk9 T lb 
m 



(4.8) 



We could ask the same question about equations of mo- 
tion when starting with the ordinary Lax pair (|4.8p satis- 
fying (|3.1|l . In this case the situation is slightly different. 
Using the explicit formula for the Christoffcl coefficients 



r/£ one can show that 



[L,M]\ = {L%, ±-p>}, 



(4.9) 



where the right-hand side is the Poisson bracket of the 
components of the Lax tensor L. Since L% = {L%, H}, 



We can thus eliminate the gamma matrices from both 
sides. Substituting the expression for the Ricci coeffi- 
cients in terms of the derivatives of components of e^, it 
is possible to show that 



1 m fm f 2l 

m I ira J 



(4.14) 



Together with p\ = {ps,, H} it leads to an equation anal- 
ogous to Eq. P~10"]) . 



{^-2^ 2 }=° 



(4.15) 



Beware however, that the frame component pa of the mo- 
mentum is not the canonical coordinate conjugated with 
x a (indeed, Pa = p n e2{x)), so (|4.15[) does not imply that 
H — 9^J> 2 is independent of the coordinates. Similarly 
to the previous case, the geodesic Hamiltonian solves this 
condition but it is an open question whether this solution 
is unique. 



9 This tensor is in some sense 'Dirac's square root' of the Lax 
tensor 114.41 1: it can be obtained from the WKB approximation of 
the Dirac equation. 
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V. LAX TENSORS AND HIDDEN 
SYMMETRIES 

Until now we have investigated rather trivial examples 
of Lax tensors, constructed from momenta and metric- 
related quantities, and hence trivially conserved along 
the geodesic motion. Such objects are present in any 
generic spacetime. Now we shall turn to examples of Lax 
tensors present in spacetimes with enhanced symmetries, 
in particular admitting hidden symmetries of Killing and 
Killing- Yano tensors. 



Lax tensor from conserved quantity 

Let us assume that, provided a given enhanced sym- 
metry of the spacetime, an additional conserved quan- 
tity E for geodesic motion is known. This, for example, 
incorporates the case of E generated from a Killing vec- 
tor £, E = £, a Pa, or E generated from a Killing tensor 
k. E = k PaPb ■ ■ ■■ Of course, L = E is a trivial one- 
dimensional Lax tensor, which can be helpful, for exam- 
ple, if one constructs a larger Lax matrix as a sum of 
independent segments. One can also simply upgrade the 
Lax tensors P~4")l or (|4.11[) as 



L E a b = Ep a p bl A E = E Pal a . 



(5.1) 



The corresponding Lax pair matrices follow from p.4[) 
or (|3.7p . If more than one such constants are known, 
we may employ the method of spectral parameters @] to 
combine the corresponding Lax tensors. 



Two kinds of Killing— Yano tensors 

The Killing-Yano (KY) tensor <p ai ...a r 1S @ an anti- 
symmetric form on the spacetime, the covariant deriva- 
tive of which is determined by its antisymmetric part, 
i.e., by its exterior derivative dip: 



V a iPa 1 ...a r = V[ a <p ai . 



(5.2) 



On the other hand, the closed conformal Killing-Yano 
(CCKY) tensor h ai ...a r is an antisymmetric form on the 
spacetime, the covariant derivative of which is deter- 
mined by its divergence £ a2 ...a r : 



D-r+ 



9a[a± £d2 ...a r 

-V h n 
1 vn"' a 2 ...a r 



(5.3) 



KY and CCKY tensors are related to each other through 
the Hodge duality: the Hodge dual of a KY form is a 
CCKY form and vice versa. 

Both KY and CCKY tensors encode the so called hid- 
den symmetries: they exist only for special metrics and 
provide a rather rich structure to the geometry, see, e.g., 
15] and references therein. 



Covariantly conserved tensors 

As one consequence of the above definitions, both KY 
and CCKY tensors define a tensorial quantity linear in 
momentum which is conserved along trajectories of the 
geodesic motion. Namely, for a KY tensor ip and a CCKY 
tensor h the quantities^ 

K = ip-p, (i = hAp, (5-4) 

are conserved along geodesies, 

l K = °> l» = °- ^ 

The velocity is given by u = ^p and the momentum p 
is conserved for the geodesic motion. Thus, thanks to 
(JOJ), the derivative of k gives ^p k (V[ k (p ai ... ar _ lt ])p l = 0. 
Similarly, Eq. (|5.3[) implies ¥rh = -^p A h which vanishes 
when wedged with another p. Of course, both equations 
in (|5.5[) are equivalent through the Hodge duality: k = *fi 
for ip = *h. 

It is interesting to observe that the same information 
encoded in the forms n and fi is also encoded in the forms 

$ = cAp, F = fi-p, (5.6) 

respectively. Indeed, employing the identity 

(a ■ p) A p + (a A p) ■ p = p 2 a , (5-7) 

which holds for any antisymmetric form a, one can re- 
construct k and fi from the quantities above 



1 



1 



2 $-p, fi=— FAp. (5.8) 
p A p l 

Since the momentum is conserved for the geodesic mo- 
tion, we thus obtained conserved quantities quadratic in 
momenta 



V V 

$ = 0, — F = 
dt dt 



(5.9) 



which, however, carry the same information as those de- 
fined in ()5.4|) . Note also the Hodge duality $ = *F for 

ip = *h. 

Finally, the form F can be rewritten using (|5.7p as 



F = (h A p) ■ p = p 2 h - (h ■ p) A p , 
which in indices reads 



(5.10) 



F. 



UhAp)-p) 

\\ i / -* / a-\_a,2...a r 



— P h j a\a'2...a r ^a±a2 ■ ■ .a r - ±nP Va r 
( l) h'a2a,3...a r nP Pa\ ■ 
( l) ha r ai ...a r -2nP Pa r -\ 



(5.11) 



P h ni 



U2...n r 1 ai 1 a2 



The dot "•" denotes the contraction, (ip ■ p) ai ...a r = </>a 1 ...a T .nJ? n - 
Since we assume automatic rising of indices, it is essentially the 
scalar product. 
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where we have introduced the projector 
P% = 6% 



P 2 p a Pb 



(5.12) 



The form F is thus (up to the prefactorp 2 ) the projection 
of the CCKY tensor h onto a subspace orthogonal to the 
momentum p. 



Killing— Yano and Lax tensors 

All the conserved tensorial quantities k, /j,, and F 
constructed from KY and CCKY tensors can be con- 
verted into Lax tensors using gamma matrices. For ex- 
ample, taking n and p, the corresponding Clifford objects 
reacFH 



Ah = l ai - ar+1 {hAp) ai ...a r+1 ■ 



(5.13) 
(5.14) 



They are covariantly conserved (|3.5|) and the correspond- 
ing Lax pair of matrices are given by (|3.7|) . As mentioned 
above, the invariants of A v span the same functional 
space as invariants of in other words, the conserved 
scalar observables generated from the Lax tensor A v are 
the same as those generated directly from and simi- 
larly for Ah and F. 

Among all conserved quantities constructed from KY 
or CCKY tensors there are special cases which do not 
need to use Clifford objects and gamma matrices. Such 
a situation occurs if the conserved tensor under consid- 
eration is of rank 2. This includes: 

(a) The CCKY tensor h a of rank 1, which is in fact a 
closed conformal Killing vector 



V /i& = t;g a b 



1 

D 



(5.15) 



The tensor fi is now of rank 2, and hence generates di- 
rectly the Lax tensor 



L a b = h a Pb - p a h b 
The Lax pair matrices are given by 

1 



— n n r a 

m 



(5.16) 



(5.17) 



cf. Eq. (|3.4[) . The only independent constant of mo- 
tion which can be obtained from this Lax tensor is 
p 2 h-P-h = p 2 h 2 - (p ■ h) 2 . 

(b) The CCKY tensor h ab of rank 2 generates the 
second-rank conserved quantity F ab given by (|5.6[) or 
(|5.1ip . The Lax tensor thus reads 



b = b b =ph b -pp tlnb - h PnPb 



(5.18) 



The Lax pair matrices are again given by (|5.17[) . In the 
special case when h ab is nondegencrate, its very existence 
guarantees complete integrability of the geodesic motion, 
see the following subsection. 

(c) The KY tensor (p abc of rank 3 generates the Lax 
tensor 



ra a n 

L b — Lp bnP 



(5.19) 



This case has been discussed in papers 0-0] • 

(d) The KY tensor ip ab of rank 2 generates the Lax 
tensor 



L a b = <S> a b = p an p nPb +p a p n p nb 



(5.20) 



The invariant generated from this Lax tensor is a function 
of the observabl J*^1p„ Vhk ab , with k a b — (p an ip bn being the 
Killing tensor of rank 2 associated with the KY tensor 

fab- 

(e) Any rank two tensor constructed from the momen- 
tum p and quantities k and p, (for various KY forms ip 
and CCKY forms h) by contractions and wedge opera- 
tion. A simple interesting example is a "partial square" 
of the quantity «, namely the Lax tensor 



L b — (f p m fbkl...nP 



(5.21) 



Trace of this Lax tensor gives the quadratic conserved 
observable p a p b k ab , where the second rank Killing tensor 
k ab is associated with the KY form <p: 



k 6 = (p tp bk i. 



(5.22) 



Another interesting possibility is to take a 'square' of 
tensor (15.61). 



La jpakl... jp 

b = -f bi 



bkl. 



(5.23) 



Depending on the rank of tp this will generate certain 
number of conserved quantities. 

To summarize, the Lax tensors build from the KY and 
CCKY forms can be very fruitful. They can generate 
plenty of functionally independent invariants and they 
are thus very useful when investigating the conserved 
quantities. 



Kerr-NUT-(A)dS spacetime 

A highly nontrivial example of the Lax tensor dis- 
cusssed above can be found in the spacetime equipped 
with a nondegenerate CCKY tensor h ab of rank 2. It 
was proved in [l8l - [20l | that the existence of such a prin- 
cipal CCKY tensor determines the form of the metric 
up to a set of metric functions of a single argument. 



11 Similar to Example II in the previous section, these Lax tensors 
can be understood as arising from the WKB approximation to 
the Dirac symmetry operators K v and studied in Il7ll . 



Here we used that p 2 is also the conserved quantity and we have 
canceled it out from tr(L 2 ). 
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These functions can be fixed by the Einstein equation 
and the resulting vacuum (with cosmological constant) 
metric describes a generally rotating black hole in arbi- 
trary number of spacetime dimensions, also called the 
Kerr-NUT-(A)dS spacetime [US]. 

In this case, the principal CCKY tensor h a b generates 
the covariantly conserved Lax tensor F = (h Ap) • p , cf. 
also (|5.18|) . which in D = 2n + e [e = 0, 1) dimensions 
gen erates n independent constants of geodesic motion 
[241427J ]. It is possible to choose these constants in such 
a way that they are quadratic in momentum and hence 
generated by rank 2 Killing tensors kfy . Namely, it was 
shown in [U| that for any parameter j3 the following iden- 
tity holds: 



p 2 det(j + ^p- 2 F) = Y,P"Po k U) P 



(5.24) 



3=0 



The left-hand side is a scalar expression, constructed just 
from the Lax tensor F% and p 2 , and hence is conserved 
for any (3. The coefficients in the /3-expansion are thus 
also conserved and can be read from the right-hand side 
of (|5.24p . They are quadratic in momentum with the km 
being Killing tensors F^l 

Moreover, the discussed spacetime admits also n + e 
explicit (Killing vector) symmetries which supply addi- 
tional n + e conserved quantities. All these D conserved 
quantities are in involution and the system is completely 
integrable. The existence of the quadratic conserved 
quantities encoded in the Lax tensor is also a start- 
ing point for showing that the Hamilton- Jacobi, Klein- 
Gordon, and Dirac equations separate in these space- 
times [ll [H HHIl. ( Sec also Sill for separability 
of certain gravitational perturbations.) 



CCKY tensor h ab , £ a = ^V^, cf. Eq. (JO). We 
will write qA a = e£ a , where q is a charge of the particle 
and e a constant combining both the charge and strength 
of the field. The motion of the charged particle is thus 
governed by the Hamiltonian: 



H = —(pa~et a )g ab (p b ~e£ b ) 



(6.1) 



The relation between the velocity and momentum can 
be read from the Hamiltonian flow (12.911: 



(6.2) 



u a = - ( P a - eC) ■ 
m 

The covariant derivative of the momenta is 
V 



dt 1 



Pa = -(y a Qip n -ec)- 



(6.3) 



Substituting (|6.2[) and using the Killing vector condition 



V a £ n = -V„£ a one gets ^p =-eK n V n £ a =-e^£ a , i.e., 



V 

dt 



(p + eO=0 



Now we can prove that the quantity 
/i = h A (p + e£) 



(6.4) 



(6.5) 



is covariantly conserved. Indeed, using (|6.4[) . (|5.3j) . and 
flSZH), we get 



V 
~dt 

It follows that 



H = (u A £) A (p + e£) = 



(6.6) 



VI. MOTION OF A CHARGED PARTICLE 



F = n-{p + eO 



(6.7) 



It was shown in |33( that the motion of a charged 
particle in the special test electromagnetic field in the 
background of Kerr-NUT spacetimcs in all dimensions 
is also completely integrable. Here we demonstrate that 
the conserved quantities quadratic in momentum can be 
found using the Lax tensor method. 

The electromagnetic field under investigation is given 
by the vector potential A proportional to the primary 
Killing vector £ of the Kerr-NUT geometrvF^ The pri- 
mary Killing vector £ is a divergence of the principal 



is also covariantly conserved. Being the tensor of rank 2, 
this is the covariant Lax tensor for our system. It is 
constructed in a similar way as the Lax tensor from the 
previous subsection, only with substitution p — > p + e£. 
The scalars generated from F can thus be read again 
from (|5.24|) . The quadratic constants of motion are 



K(j) = (p a + e£ a )(Pb + e£ b )k 



ab 



(6.8) 



These differ from the constants Ku\ introduced in [33(, 
but only by terms linear in momentum which are 
also conserved: 



13 For j = the Killing tensor reduces to the metric, = g ab . 
The n-th Killing tensor vanishes in even dimensions, whereas it 
is reducible to a square of one of the Killing vectors in odd di- 
mensions. For j = 0, , n — 1 the Killing tensors are irreducible, 

giving thus n quadratic constants of motion. 

14 The condition of vanishing electric current for such an electro- 
magnetic field requires the cosmological constant to be set equal 
zero, cf. 1331. 



K U) = (Po ~ e £a)(Pfo - eib)kf = K {j) - 4e L {j) , 



(6.9) 



The conservation of L(j\ follows from the fact that 
^U) = ^{f)^ n are billing vectors [U]. It was demonstrated 
in [33l ] that the conserved quantities and L(j\ are all 
in involution. 
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VII. SUMMARY 



Appendix A: Derivatives on the phase space with 
cotangent bundle structure 



The Lax pair formalism provides an elegant and ef- 
fective description of special dynamical systems with en- 
hanced symmetries. In particular, the existence of the 
Lax pair, defined by Eq. (|1.2[) . enables one to generate 
constants of motion by simple algebraic operations, e.g., 
Eq. STfy . 

In this paper we have provided an alternative, covari- 
ant formulation of the Lax formalism. This is based on 
the covariant (Clifford) Lax tensor, where the Lax equa- 
tion is formulated as a covariant conservation of this ten- 
sor, Eqs. (|3.2[) and (|3.5j) . In both instances the existence 
of the Lax tensor enables one to generate constants of 
motion, which are determined as invariants constructed 
from the object, e.g., Eq. (|1.5p . We have further demon- 
strated that the ordinary Lax pair matrices follow from 
the covariant Lax formalism, the relation being given by 
Eqs. ([53J and (|577)l . 

To illustrate the derived formulas, we have concen- 
trated on the problem of particle motion in curved space- 
time. In this case we were able to provide a number of 
examples of (Clifford) Lax tensors. In particular, we have 
concentrated on manifolds with enhanced symmetry, ad- 
mitting hidden symmetries of Killing- Yano tensors, in 
which case the examples of Lax tensors are highly non- 
trivial. One of the Lax tensors discussed was proved re- 
sponsible for complete integrability of geodesic motion 
in rotating black hole spacctimcs in all dimensions just a 
few years ago. We have demonstrated for the first time 
that the conserved quantities for motion of a charged par- 
ticle in the aligned test electromagnetic field on the same 
(vacuum) black hole background can be also generated 
using the Lax tensor (|6.7p . 

It remains an interesting open question whether any of 
the Lax tensors discussed here will find further physical 
applications in the future. 
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In this Appendix we discuss the structure of the cotan- 
gent bundle phase spacc0 P = T*M in more details. 
We show that the covariant derivative on the configura- 
tion space M induces a covariant splitting of the phase- 
space quantities into quantities related to the configu- 
ration space. It naturally replaces standard coordinate 
expressions in a coordinate independent way. 



Derivatives along position and momentum directions 

First, we introduce covariant partial derivatives of a 
scalar observable along position and momentum direc- 
tions. The derivative in momentum direction f a (i.e., 
changes along a curve p a — > p a + sf a , x fixed) is simple, 
since the space of momenta at fixed x is linear. We define 



fa 



dF _ d 
dp a de 



F(x,p + ef) 



e=0 



(Al) 



Thanks to ultralocality in f a we can tear off f a to 
obtain derivative operator (with one contravariant 
configuration-space index) acting on scalar phase-space 
observablcs. Such a derivative operator also defines a 
mixed tensor J^- e TP ® TM. It is actually the ten- 
sor identifying the tangent space of the cotangent fibre 
T(T*M) with the cotangent fibre T*M itself. 

The derivative along a position direction with momen- 
tum fixed is more involved. Moving from one position 
to another one changes the cotangent fibre and it has 
to be clarified what "fixed momentum" means. A nat- 
ural solution is given in terms of a spacetime covariant 
derivative^ V. The covariant derivative defines "fixed 
momentum" to be the parallel-transported momentum. 
Let x £ be a spacetime curve in the u a direction and p e 
be the parallelly transported momentum along this curve. 
Then we can write 



a VaF d _ 

u -5— = —F(x E ,p E 
ox de 



e=0 



(A2) 



Again, it defines the phase space derivative ^ with one 
covariant configuration-space index and the mixed tensor 
v A 

-rg: G TP ® T*M. The last quantity is the tensor which 
makes a horizontal lift of the configuration-space vector 

u a to the horizontal phase-space vector u a -g^ , see Fig. [1] 



15 We use capital Latin indices for the phase-space tensors and in 
this Appendix we write these indices explicitly. The material 
presented here partially follows and partially generalizes the Ap- 
pendix of [2(1 , 

16 It can be an arbitrary covariant derivative. Of course, in most 
cases it is useful to chose the metric derivative. 
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The action of these derivatives on an observable of the The completeness relation reads 



form A(x,p) = a ' (x)p a Pb • • • is 



dx 
dA 

dp n 



V n a ab -{x)p a p b . 
a nb -(x)p b --- + 



'{x)Pa ■ 



(A3) 
(A4) 



The action on a general observable can be written ex- 
plicitly using linearity, the Leibnitz product rule, and 
the chain rule. 

The mixed tensors and -g— are covariant general- 
ization of the phase-space coordinate vector^ -§^s an d 
jj^- associated with the canonical coordinates (x a , p a ). 
Therefore, one can expect that natural symplectic quan- 
tities can be written using these covariant tensors. 
Namely, the inverse symplectic structure Q~ 1AB (such 
that ^Ian^ 1BN = the Poisson brackets, and the 
Hamiltonian flow Xjj = il~ 1AN d^H are: 

dx dp n dp n dx " kl dpk dpi 
sprn V n FdG OF V„G dF dG 

H dp n dx 



-IAB 



dx dp n " dpk kl dpi 



(A7) 



Here, T$ is the torsion of V. For vanishing torsion T = 
the expressions resemble the standard coordinate formu- 
lae. 



D n A x 



dx 



^APn 



dp n 



'A ■ 



(A9) 



The symplectic structure Qab and the symplectic poten- 
tial 9a = —p n dAX n have a familiar form (except for the 
anomalous torsion term): 

Vab = Dlx V B p n -V A p n D£x-p n T£ D\x D l B x , (A10) 
e A = -p n D n A x. (All) 

These dual mixed tensors define also the splitting of a 
phase space vector X A into configuration-space quanti- 
ties u a and f a discussed in Sec. [TTI near Eq. (|2.5[) . Such 
splitting is based on the observation that the covariant 
derivative splits the the tangent fibre T^piP into hori- 
zontal and vertical subspacesEl The configuration direc- 
tion u a is the projection of X A onto spacetime (encoding 
thus the horizontal part) and the momentum direction 
f a is the vertical part of X A . These two parts can be 
written using Dx and Vp, see Fig. [TJ 

Indeed, the tensor Dx is just the differential of the 
projection map x : [x,p] — > x from the phase space to the 
configuration space. Given a phase-space vector X A , its 
shadow on the configuration space thus is 



= X J 



D a A x 



(A12) 



On other hand, the quantity Vp is the projector of a 
phase space vector X A on its vertical part fa * 



Dual quantities and splitting of the phase-space 
direction 



We can also construct dual quantities D^ x and V4 Pn 
which corresponds to coordinate forms d,AX a and d A p a - 
We require the duality conditions 



dx ° nX 



qN 

a VNPb 

dp a 



°b ) 



dx 



y NPb = , ^— D b N x = 



QN 

dpa 



(A8) 



fa = X V APa 



(A13) 



The completeness relation (|A9[) together with (|A12I) 
and (|A13|1 gives 



X' 



V A d A 



dx 



dp n 



(A14) 



cf. Fig. [TJ and the derivative of scalar the observable 
F(x,p) along the phase-space direction X A is thus 



vA j n ^nF dF 

X A d A F = u n —— + /„ — 
dx dp n 



(A15) 



17 The notation here is a bit mistreating the difference between co- 
variant nature of the mixed tensors and coordinate-dependent 
nature of the coordinate tensors. For example, both indices in 

V A 

-a 3 - are tensor indices and could be understood as abstract in- 

OX 

dices, independent of the chosen coordinates. Similarly for the 

r) A 

phase-space index A in the coordinate vector ^h - - However, the 
index n here is not a tensor index, it just labels which coordi- 
nate tensor we are choosing. To make this distinction clear, one 
should distinguish the abstract and coordinate indices as, e.g., in 
[34j . We decided not do so and let the reader distinguish tenso- 
rial and coordinate indices based on the context. In general, in 
expressions not involving explicitly chosen coordinates all indices 
are tensorial. If the coordinates are involved, the combinations 
x a and p a usually indicate the coordinate indices. 



The covariant derivative induced on the phase space 

The covariant derivative of the phase space fields with 
configuration space indices introduced in Sec. UU can be 
understood as a generalization of the formula (|A15[) to 



The horizontal subspace of Ty xp ^P gives the directions of 
parallel-transported momenta, the vertical subspace is tangent 
to the fibre TIM. 
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P = T'M 



on tensor field 




The action of the derivatives 5* and 

OX 

A%---(x,p) = a%' PkPi ■ ■ ■ is analogous to rules (|A3j) 
and (|A4|) . just with additional indices involved. 

The covariant derivative along a general phase-space 
direction X A (split as in (|A14[0 acting on a field 
A%'"(x,p) is then given by generalization of (|Al5j) : 



° B S^» + ^^--- (Al7) 



This we have already mentioned in Eq. (|2.7p . In par- 
ticular, the derivatives of the pure configurations field 
A^'"{x) and of the momentum field p a are given by rules 
(i) and (ii) in Sec. HI 

Since the dependence on X is ultralocal, it is possible 
to define also the covariant differential V/vA^- ', namely 



V N A a bZ 



ryn„ YlL A a — 
° nX dx A " 



H... 



^NPn ^-^6 " 
Op n " 



(A18) 



Relation to canonical coordinates 



u=X-Dx 
M 



FIG. 1. Splitting of the phase-space direction X. In this 
figure we illustrate the splitting of the phase-space direction 
X into its position and momentum parts. Let [a; £ ,p c ] be a 
curve starting at [x,p] to which X is tangent. Its correspond- 
ing position and momentum parts are x £ and p £ . Dashed 
curve p e is the parallel transport of the initial momentum p 
along x e - Short-dashed curve p £ laying in the cotangent fi- 
bre TJM is obtained as the parallel transport of p e along x e 
back to the point x. Various vectors defined in the text are 
tangent vectors to these curves. The phase-space vector X is 
tangent to [x £ ,p £ ], its position direction u is tangent to x e , 
the momentum part / is the derivative of p e understood as a 
cotangent vector. In other words, / is the covariant derivative 
of p £ along x e , i.e., along the direction u. The formula (|A14|) 
splits X into u ■ which is tangent to p e , and into / • J^, 



which is tangent to p s 



tensor fields. The derivative acting in the momen- 
tum directions has the same definition (|A1[) . The defi- 
nition of the derivative 5 s - in the configuration direction 



changes just employing the covariant derivativd I in the 
definition 



y*K.:: 

dx 



—A%;;;(x e ,p e ) 



e=0 



(A16) 



19 Since A^'" (x e , p e ) is just e-dependent, it can be understood as 
a tensor field along the configuration curve x e and just the stan- 
dard covariant derivative in the configuration space is involved 
on the right-hand side of the definition I IA16I I. 



In the definitions above we have not used any specific 
choice of the coordinates. However, the introduced for- 
malism can be easily accommodated to such a choice. 
If we choose configuration-space coordinates x a , one can 
define the "coordinate derivative" d by conditions 



ddx b = , d 



d 
dx b 







(A19) 



It is a torsion-free covariant derivative (of course, de- 
pending on the choice of coordinates). The difference 
tensor between V and d is given by the connection coef- 
ficients r™ b . 

We can use the coordinate derivative d instead of V 
in all expressions above. It leads to standard coordinate 
expressions: the equation (jAlOj) reduces to (|2.1|) . the ex- 
pression (|A6|) to (|2.2I) . etc. In the coordinate case we 
also use the more common notation 7A7 instead of Jr-. 

ax a ox 

The induced coordinate derivative on phase space along 
a phase-space direction X A is denoted just by dot 



d 

A a... _ ±_ Aa 

b - ~ de b 



(A20) 



Splitting of the phase-space direction X A depends on 
the choice of the covariant derivative. For the coordinate 
derivative this splitting gives: 



u a = £ A D%x 



Pa = £ 8 A Pa 



(A21) 



Clearly, p a are just derivatives of components of p along 
X A - which justifies the dot notation. 

Let [x e ,p e ] be a phase-space curve with tangent vec- 
tor X. Then, f a = X A V A Pa can be understood as the 
standard configuration-space covariant derivative of p e 
along the spacetime curve x s , cf. Fig. [TJ Similarly, 
p a = X A d A Pa is the coordinate derivative of p e along x e . 
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We can thus use the ordinary relation between covariant 
and coordinate derivatives to obtain 



fa =Va-U k Tl aP i , 



(A22) 



cf. (|2.6I) . Here, the velocity u k is tangent to the curve x e . 
Tearing off the phase-space vector X A , we obtain 

V APn = d APn - D\x Yl n pi . (A23) 

The duality relations (|A8|) imply 



• n _^ I j^k _^ 

dx ~ dx n Pk nl dpi 
The action on a scalar observable reads 



dF 



(A24) 



(A25) 



<9x dx n 

The generalization to the action on tensor fields adds 



only standard terms for each tensor index: 

g A a.. 



v, 



n Aa... _ ^ Jl b... , -pk u ^b.. 

' dx" +PkKl d Pl 



rltl A k. . . "pifc ylCI... 



(A26) 



rib k... 



If we substitute ()A22j) and (|A26|) into (|A17|) . we obtain 
the coordinate expression for the covariant derivative on 
the phase space 



V • i 

— A a ■■■ — A a - - 4- n n r" A k - 



Here, the coordinate derivative A%- - splits as 



Mr = u n 0At:: 



dA a... 



dx r - 



dpn 



(A27) 



(A28) 



P. Lax, Comm. Pure Appl. Math 21, 467 (1968) 
O. Babelon, D. Bernard, and M. Talon, Introduction to 
classical intcgrablc systems (Cambridge University Press, 
Cambridge, 2003) 

K. Rosquist(1994), |arXiv:gr-qc/9410011 [gr-qc]| 

K. Ros quist and M. Goliath, IGen.Rel.Grav.l 30. 1521 

(1998) , |arXiv:gr-qc/9707003 [gr-qc]| 

M. Ka rlovini and K. Rosquist, IGen.Rel.Grav.l 31, 1271 

(1999) , |arXiv:gr-qc/9807051 [gr-qc]| 

D. Bale anu and A. K. Karasu, Mod.P hys.Lett. A14, 2597 

(1999) , [arXiv:gr-qc/0004024 [gr-qcfl 

D. Bal eanu and S. Baskal, Mod.P hys.Lett. 15, 1503 

(2000) , |arXiv:gr-qc/0104011 [gr-qc]| 
K. Yano, Ann. Math. 5 5, 328 (1952) 
U. Semmelmann(200 2), |math/0206117] 

G. Papadopoulos, Class. Q uant. Grav.| 25, 105016 
(2008) , |arXiv:0712.0542 [hep-thJT 

V. P. Frolov and D. Kubiznak, Class. Quantum Grav. 25, 
154005 (2008), [arXiv:0802. 0322 [hep-th]| 
Y. Yas ui and T. Houri,|Prog.Theor.Phys.Suppl.| 189, 126 
(2011), [arXiv:1104.0852 [hep-th]| 

O. Santillan, |J.Math Phy'sTT 53, 043509 (2012), 



|arXiv:1108.0149 [hep-th][ 
M. Visinescu and G. E7 Vilcu, SIGMA 8, 058 (2012), 

|arXiv:1205.6036 [hep-th]] 

M. Cariglia(2012), |arXiv:1209.6406 [gr-qc]| 
M. Cariglia, P. Krtous, and D. Kubiznak, 
D84, 024004 (2011), |arXiv:1102.4501 [hep-th] " 
M. Cariglia, P. Krtous, and D. Kubiznak, 
024008 (2011), [arXiv: 1104.4123 [hep-th]' 



Phys.RevT] 



D84, 



T. Hou ri, T. Oota, and Y. Yasui, |Phys. Lett. B656, 214 



(2007) , |arXiv:0708.1368 [hep-th]] 
P. Krtous, V. P. Frolov, and D. Kubiznak, Phys. Rev. 



Phys.Rev. 



[20 
[21 
[22 

[23 
[24 
[25 
[20 

[27 
[28 
[29 
[30 
[31 
[32 
[33 
[34 



D78, 064022 (2008), [arXTv:0804.4705 [hep-th]| 

T. Houri, T. Oota, and Y. Yasui, Phys. Lett. B666, 391 



(2008), arXiv:0805.0838 [hep-th] 



172, 



R. C. Myers and M. J. Perry, Ann. Phys. (N.Y. 
304 (1986) 

G. W. Gibb ons, H. Lii, D. N. Page, and 

C. N. Pope, |Phys Rev. Lett.| 93, 171102 (2004), 
|arXiv:hep-th/0409155| 

W. Chen, H. Lii, and C. N. Pope, Class. Quantum Grav. 
23, 5323 (2006), [arXiv:hep-th/0604125| 

D. N. Page, D. Kubiznak, M. Vasudevan, and P. Krtous, 
Phys. Rev. Lett. 98, 061102 (2007), [hep^th/0611083| 

P. Krtous, D. Kubiznak, D. N. Page, and V. P. Frolov, J. 



High Energy Phys. 02, 004 (2007) , [ai^v:hep-th/0612029 
P. Krtous, D. Kubiznak, D. N. Page, and 
M. Vasudevan, |Phys. RevJ D76, 084034 (2007), 

|arXiv:0707.0001 [hep-th]T 

T. Hou ri, T. Oota, and Y. Yasui , [J. Phys.| A41, 025204 



(2008) , |arXiv:0707.4039 [hep-th 



V. P. Frolov, P. Krtous, and D. Kubiznak, J. High Energy 
Phys. 02, 005 (2007), [hep^"th/0611245| 
A. Ser gyeyev and P. Krtous, |Phys. Rev] D77, 044033 
(2008) , |arXiv:071 1.4623 [hep-thX~ 

T. Ootaaiid Y. Yasui, [Phys. Lett.] B659, 688 (2008), 
arXiv:0711.0078 [hep-th][ ~ 

H. K. Kunduri, J. Lucietti, and H. S. Reall, Phys. Rev. 
D74, 084021 (2006), |axXTv:hep-th/0606076| 
Oota and Y. Yasui, |Int.J.Mod.Phys? 



T 



A25, 3055 



(2010), |arXiv:0812. 1623 [hep-th]| 
V. P. Frolov and P. Krtous, |Phys. Rev. DlD83, 024016 



(2011),|arXiv:1010.2266 [hep-th 



R. Penrose and W. Rindler, Spmors and Space-Time 
(Cambridge University Press, Cambridge, England, 1984, 
1986) 



